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(54) Method and apparatus for three-dimensional image rendering of body organs 



(57) A method for mapping a structure in a body of 
a subject includes capturing a three-dimensional (3D) 
image of the structure comprising diagnostic informa- 
tion, and geoeraUng a 3D geometrical map of the struc- 
ture using a probe inserted into the structure. The image 



FIG. 1 



Is registered with the map, such that each of a plurality 
of image points in the Image is Identified with a corre- 
sponding map point in the map. The map is displayed 
such that the diagnostic information associated with 
each of the Image points is displayed at the correspond- 
ing map point. 
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Description 

FIELD OF THE INVENTION 

s {0001 J The present invention relates generally to systems and methods for three-dimensional mapping and recon- 
struction, and specifically to mapping and reconstruction of the interior of body organs, such as the heart. 

BACKGROUND OF THE INVENTION 

10 [0002] Various methods of diagnostic imaging are known in the art. Methods used for imaging the heart, for example, 
include fluoroscopy, angiography, echocardiography, computed tomography (CT), magnetic resonance imaging (MRl). 
positron emission tomography (PET) and single photon emission tomography (SPECT). Many of these methods pro- 
duce three-dimensional (3D) Image Information, which can then be rendered for viewing in the form of parallel slices 
through the heart, or as a pseudo-3D display on a video monitor. In order to administer treatment, the treating physician 

J5 must build a 3D picture in his or her mind based on the two-dimensional pictures that are displayed. The transposition 
Is particularly trtclcy when therapy is to be administered inside the heart, such as local electrical ablation of aberrant 
electrical pathways, or laser myocardial revascularization. 

[0003] It Is also known In the an to map the heart using a mapping probe, typically a catheter, Inside the heart 
chambers. Exemplary methods and devices for this purpose are described in U.S. Patents 5.471,982 and 5,391 ,199 
.20 and in PCT patent publications WO94/06349, WO96/05768 and WQ97/249B1, whose disclosures are incorporated 
herein by reference. U.S. Patent 5,391 ,199, for example, describes a catheter that includes both electrodes for sensing 
cardiac electrical activity and miniature coils for determining the position of the catheter relative to an externally-applied 
magnetic field. Using this catheter a cardiologist can collect data from a set of sampled points In the heart within a 
short period of time, by measuring the electrical activity at a plurality of locations and determining the spatial coordinates 
55 of the locations. Locations of the mapping catheter within the heart can be superimposed on a 3D reconstruction of an 
Image of the heart, such as an ultrasound Image, acquired prior to or during the catheter study. Color codes are used 
to represent electrical activity sensed by the catheter. 

[0004] U.S. Patent 5,738.096, whose disclosure is incorporated herein by reference, describes methods for geomet- 
rical mapping of the endocardium based on bringing a probe into contact with multiple locations on a wail of the heart. 

so and determining position coordinates of the probe at each of the locations. The position coordinates are combined to 
form a map of at least a portion of the heart. Once the position of the catheter is known, external sensors can be used 
to provide local physiological values of heart tissue adjacent to the tip of the catheter. For example, If the catheter 
incorporates a radioactive marker suitable for SPECT, local functional information can be gleaned from a SPECT image. 
Yet another example is determining local perfusion from Doppler-ultrasound images of the coronaries, from nuclear 

3S medicine images or from X-ray or CT angiography, and overlaying the perfusion map on the geometrical map. The 
image of the catheter In the perfusion map can be used to align the perfusion map and the geometrical map. Alterna- 
tively, the alignment may be carried out using fiducial marks or anatomical reference locations, either automatically or 
manually. 

[0005] Further methods for creating a three-dimensional map of the heart based on these data are disclosed, for 
40 example, In European patent application EP 0 974 936 and in EP-A-0 974 936 whose disclosure Is incorporated herein 
by reference. As indicated In these applications, position coordinates (and optionally electrical activity, as well) are 
initially measured at about 10 to 20 points on the interior surface of the heart These data points are generally sufficient 
to generate a preliminary reconstruction or map of the cardiac surface to a satisfactory quality. The preliminary map 
is preferably combined with data taken at additional points In order to generate a more comprehensive map. 

SUMMARY OF THE INVENTION 

[0006] It is an object of some aspects of the present invention to provide improved methods and apparatus for map- 
ping and visualization of internal body structures, and particularly of the heart. 
so [0007] It is a further object of some aspects ol the present invention to provide improved methods and apparatus for 
administering local treatment of pathological conditions within Ihe heart. 

[0008] In preferred embodiments of the present invention, a position-sensing catheter Is used to generate a 3D 
geometrical map of the internal surface of a heart chamber of a subject. A 3D diagnostic image of the heart is captured 
in conjunction with generating the 3D map, typically either before or concurrently with the mapping. The image and 
35 map are brought into mutual registration, and diagnostic information from the image, such as perfusion information, is 
then marked on the 3D map, preferably in the form of color coding. Based on the combined diagnostic and geometrical 
information, a physician operating the catheter is able to identify and visualize areas of the heart that are in need of 
treatment, due to low perfusion, for example. The physician preferably uses the catheter to apply a local invasive 
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therapy, such as laser revascularization, to specific points that are located using the color-coded 3D map. Alternatively, 
a local diagnostic technique, such as a biopsy, may be performed at such specific points. 

[00091 There Is therefore provided, in accordance with a preferred embodiment of the present Invention, a method 
for mapping a structure in a body ot a subject, Including: 

5 

capturing a three-dimensional (3D) image of the structure including diagnostic information; 
generating a 3D geometrical map of the structure using a probe Inserted Into the structure; 
registering the image with the map, such that each of a plurality of image points in the image Is identified with a 
corresponding map point In the map; and 
w displaying the map, such that the diagnostic information associated with each of the image points is displayed at 
the corresponding map point. 

.,; . [0010] In a preferred embodiment, the diagnostic information is related to blood flow in the structure, wherein the 
diagnostic Information includes local perfusion data. In other preferred embodiments, the diagnostic information In- 

J5 eludes metabolic data, or Is related to uptake of a substance in tissue of the structure, or is related to motion of the 
structure. 

[0011 J Preferably, generating the geometrical map Includes bringing the. probe Into contact wtth the structure at a 
multiplicity of locations on the structure, and recording position coordinates of the probe at the locations, wherein 
recording the position coordinates includes determining the coordinates using a position sensor in the probe. 

20 [001 2] Preferably, registering the image with the map includes applying a transformation to at least one of the image* 
and the map so that following the transformation, the image and the map have a common axis and a common scale. 
Further preferably, registering the image with the map includes dividing the Image into a plurality of parallel planar 
slices, perpendicular to the axis and mutually spaced along the axis, wherein the plurality of image points are located 
in the slices. More preferably, registering the image with the map includes finding an axial coordinate of each of the 

25 slices and an angular coordinate of each of the image points located in each of the slices, and Identifying each of the 
image points with the map point having the same axial and angular coordinates. Most preferably, the structure includes 
a wall defining a cavity, and identifying each of the image points with the map point includes finding, at the axial and 
the angular coordinate, the image point that is within a section of the wall. 

[001 3] Preferably, displaying the map Includes coloring the map to reflect the diagnostic information. 

30 [0014] In a preferred embodiment, the method Includes performing a medical procedure on the structure guided by 
the diagnostic information displayed on the map. Preferably, performing the medical procedure includes using the 
probe to perform the procedure locally at locations selected on the geometrical map, and the method includes marking 
on the geometrical map the locations at which the procedure was performed. Additionally or alternatively, performing 
the medical procedure includes performing a therapeutic procedure, wherein the diagnostic information relates to local 

35 blood flow in the structure, and wherein performing the therapeutic procedure includes performing a procedure for 
improving the local blood flow. Alternatively, performing the medical procedure includes performing a diagnostic pro- 
cedure. 

[001 5] Preferably, the structure includes a heart of the subject, and generating the geometrical map includes mapping 
an endocardial surface in a ventricle of the heart. 
40 [0016] There is also provided, In accordance with a preferred embodiment of the present invention, apparatus for 
mapping a structure in a body of a subject, including: 

an Imaging device, adapted to capture a three-dimensional (3D) image of the structure including diagnostic Infor- 
mation; 

45 a probe, adapted to be Inserted into the structure, so as to generate a 3D geometrical map of the structure; 

a processor, coupled to the probe and to the Imaging device, and adapted to register the Image with the map, such 
that each of a plurality of image points in the image is identified with a corresponding map point In the map; and 
a display, coupled to be driven by the processorto display the map, such that the diagnostic information associated 
with each of the image points is displayed at the corresponding map point. 

so 

[0017] The present invention will be more fully understood from the following detailed description of the preferred 
embodiments thereof, taken together with the drawings in which: 

BRIEF DESCRIPTION OF THE DRAWINGS 

55 

[0018] 

Fig. 1 is a schematic, pictorial Illustration of a system for imaging, mapping and treatment of the heart, in accordance 
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with a preferred embodiment of the present invention; 

Fig. 2 Is a flow chart that schematically illustrates a method for Imaging, mapping and treating the heart, In accord- 
ance with a preferred embodiment of the present invention; 

Rg. 3 is a schematic representation of a map of a chamber of the heart, in accordance with a preferred embodiment 
of the present invention; 

Rg. 4 is a simplified geometrical representation o! the map of Fig. 3, showing coordinates used in registering the 
w map with an image of the heart, in accordance with a preferred embodiment of the present invention; 

Rg. 5 is a schematic, exploded view of a 3D image of the heart, represented as a stack of parallel slices through 
the heart, In accordance with a preferred embodiment of the present invention; 

is Rg. 6 shows the slices of Fig. 5 arrayed slde-by-side, illustrating registration of the slices with the 30 map of Rg. 

3, in accordance with a preferred embodiment of the present Invention; and 

Rg. 7 Is a schematic representation of the map of Fig. 3, after coloring of the map with diagnostic Information from 
the Image of Figs. 5 and 6, in accordance with a preferred embodiment of the present invention. 

20 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

[001 9] Fig. 1 is a schematic, pictorial illustration of a system 20 for three-dimensional geometrical mapping, imaging 
and treatment of e heart 24 of a subject 26, in eccordance with a preferred embodiment of the present Invention . System 
25 20 comprises an elongate probe, preferably a catheter 30, which Is inserted by a user 22 through a vein or artery of 
the subject into a chamber of the heart 

[0020J Catheter 30 preferably comprises at least one position sensor (not shown in the figures), most preferably 
located near the catheter's distal tip. The position sensor preferably comprises an electromagnetic sensor, which is 
mounted within ihe catheter by any suitable method, tor example, using potyurethane glue or the like. The sensor is 

30 electrically connected to an electromagnetic sensor cable, which extends through the catheter body and into a control 
handle of the catheter. In the control handle, the wires of the sensor cable are connected to a circuit board (not shown), 
which amplifies the signals received from the electromagnetic sensor and transmits them to a computer housed in a 
console 34 . in a form understandable to the computer. Because the catheter is designed for single use only, the circuit 
board preferably contains an EPROM chip, which shuts down the circuit board after the catheter has been used. This 

35 prevents the catheter, or at least the electromagnetic sensor, from being used twice. 

[00211 To use Ihe electromagnetic sensor, subject 26 is placed in a magnetic field generated, for example, by situating 
under the patient a pad containing field generator coils 2B for generating a magnetic field, driven by driver circuits 32. 
A reference electromagnetic sensor (not shown) Is preferably fixed relative to the patient, e.g., taped to the patient's 
back, and catheter 30 containing its sensor is advanced into heart 24. The sensor preferably comprises three small 

to coils,' which in the magnetic field generate weak electrical signals Indicative of their position In the magnetic field. 
Signals generated by both the fixed reference sensor and by the sensor in the heart are amplified and transmitted to 
console 34, which analyzes the signals and then displays theresults on a monitor 36. By this method, the precise — r 
location of the sensor in the catheter relative to the reference sensor can be ascertained and visually displayed. The J 
sensors can also detect displacement of the catheter that is caused by contraction of the heart muscle. 

45 [0022] Suitable electromagnetic sensors for the purposes of the present Invention are described, for example, in the 
above-mentioned U.S. Patent 5,391 ,1 99 and PCT patent publication WO 96/05768. A preferred electromagnetic map- 
ping sensor is manufactured by Biosense Ltd. (Tirat Hacarmel, Israel) and marketed under the trade designation NOGA. 
Some of the mapping features of catheter 30 and system 20 are implemented in the NOGA-STAR catheter marketed 
by Biosense Webster, Inc., and in the Biosense-NOGA system, also marketed by Biosense Webster, Inc. Further 

so aspects of the design of catheter 30 and ol system 20 generally are described In EP-A-1 125 549 which is Incorporated 
herein by reference. Using such sensors, system 20 achieves continuous generation of six dimensions of position and 
orientation information with respect to catheter 30. Alternatively, the sensors used in catheter 20 may comprise other 
types of position and/or coordinate sensors, as described, for example, in U.S. Patent 6,391,199, 5,443,489 or 
5,515,853, or in PCT publication WO 94/04938 or WO 99/05971, or substantially any other suitable type of position/ 

55 coordinate sensing device known in the art. 

[0023] As noted above, catheter 30 is coupled to console 34, which enables the user to observe and regulate the 
functions of the catheter. Console 34 Includes a processor, preferably a computer with appropriate signal processing 
circuits (which are typically contained Inside a housing of the computer). The processor is coupled to drive display 36. 
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User 22 brings the distal tip of catheter 30 into contact with multiple points on the endocardial surface of heart 24. and 

the position coordinates are recorded at each point. The information derived from this analysis is used to reconstruct J 
a three-dimensional geometrical map 38 of the endocardial surface of heart 24. 

[0024] System 20 also comprises a diagnostic imaging unit 48, such as an echo Doppler unit. SPECT, PET, MRI, 

5 CT or other imaging unit known in the art. Unit 48 is used to capture a 3D diagnostic image of heart 24. preferably 
while user 22 is mapping the heart using catheter 30. Alternatively, the diagnostic image Is captured before beginning 
the mapping, and unit 48 may, In this case, be separate from the other elements of system 20. Diagnostic data from 
the image captured by unit 48 are superimposed on map 38, using methods described herelnbelow. Depending on the 
type and configuration of unit 48, a wide range of different diagnostic data may be represented In the Image, such as J 

10 perfusion, metabolic factors, uptake of markers, heart wall motion or thickness, and/or other anatomical or electrical J 
properties, as are known in the art. The image can also be limed to represent different phases in the heart cycle. ^ 
[0025] Typically, system 20 Includes other elements, some of which are not shown in the figures for the sake of 
„ simplicity. In the present embodiment, the system preferably includes a laser console 49, which is used in performing 
direct myocardial revascularization, as described, for example, in PCT patent application PCT/IL97/00O11 and in U.S. 

is patent application 09/1 09,820, which Is assigned to the assignee of the present patent application and whose disclosure 
Is incorporated herein by reference. Console 49 Injects laser energy Into a suitable waveguide (not shown) within 
catheter 30. The waveguide conveys the energy to the distal tip of the catheter, where It Is applied to revascularlze 
areas of the myocardium suffering from low perfusion. Alternatively, the system may Include other therapeutic elements, 
as are known in the art, particularly elemenls for delivering local treatment In the heart, such as a radio- frequency 

20 driver coupled to an ablation electrode on catheter 30; an ultrasound generator coupled to high-power transducer in 
the catheter, for ultrasonic ablation of the endocardium; or a supply of a therapeutic agent, such as growth factors for 
angiogenesis, coupled to an Injection needle in the catheter. Still further alternatively, the system may include invasive 
diagnostic elements, such as biopsy forceps that are operated through catheter 30. ^ — ' 
[0026] Other elements that may be comprised in system 20 are described, for example, In EP-A-0 974 936. >J { \S 

25 [0027] Typically, system 20 includes an ECG monitor (not 6hown), coupled to receive signals from one or more body \ ^ ^ [> 
surface electrodes, so as to provide an ECG synchronization signal to console 34. As mentioned above, the system X 
preferably also includes a reference position sensor, either on an externally-applied reference patch attached to the 
exterior of the patient's body, or on an internally-placed catheter, which is inserted into heart 24 and maintained In a 
fixed position relative to the heart. By comparing the position of catheter 30 to that of the reference catheter, the 

30 coordinates of catheter 30 are accurately determined relative to the heart, Irrespective of heart motion. Alternatively, 
any other suitable method may be used to compensate for heart motion. 

[0028] Fig. 2 is a flow chart that schematically illustrates a method for imaging, mapping and treatment of heart 24 
using system 20, in accordance with a preferred embodiment of the present invention. At an imaging step 50. a diag- 
nostic image of heart 24, such as a SPECT image, is captured. Preferably, although not necessarily, the image is 

35 captured while catheter 30 is already located inside the heart. The catheter is used to generate geometrical map 38, 
at a mapping step 52. Suitable mapping techniques forthis purpose are described In the above-mentioned U.S. Patent 
5,738,096. The above-mentioned European patent application EP 0 974 936 describe accurate methods for creating 
the map itself based on the data gathered using catheter 30. The image captured at step 50 and the map created at 
step 52 are then registered one with the other, at a registration step 54. 

40 [0029] Figs. 3 and 4 are schematic representations of map 38 generated by system 20 at step 52, illustrating a 

method used at registration step 54, in accordance with a preferred embodiment of the present invention? Fig. 3 is a /\. 0\ JJ- 
wire frame rendition of the map, representing the left ventricle of heart 24. For the purposes of step 54. a longitudinal - "!/ 
axis 72 is drawn through the map, passing through an apex 74 of the ventricle. Preferably, the axis and apex are found '1 
automatically by console 34. Alternatively or additionally, these or other features of the map are Identified manually byj 

45 user 22. Y 
[0030] Fig. 4 Is a simplified geometrical representation of a surface 80 of map 38, generated for the purpose of 
registration with a diagnostic Image of heart 24. Surface 80 corresponds to an approximate locus of the endocardium] m 
of the heart, as determined from map 38. A coordinate system is defined in which each point 82 on surface 80 Is 
represented by a distance R from apex 74 and an angle a relative to a downward direction 84 (i.e., the direction pointing 

50 toward the feet of subject 26). A L y^J 

[0031] In order to register the diagnostic image with map 38, axis 72 and apex 74 are Identified In the Image, as well, r\V\c+\->™ 
and are aligned with the axis and apex of the map. The identification is preferably automatic but may, alternatively or f" JT 

additionally, be carried out or assisted by user 22. Other landmarks and/or fiducial marks In the heart can also be used j^-* vj I tA 
in performing the alignment. The scale of the image is adjusted so that its dimensions match those of the map asj 

55 closely as possible. For many types of diagnostic images, such as perfusion maps, the resolution of the diagnostic 
Information is low, so that imprecision of as much as 10 mm in mutual registration can be tolerated. When higher 
resolution is required, the registration of the diagnostic Image with the geometrical map may be improved using methods 
of automatic registration such as those described in Appendix A. These methods are optional and are not essential to 
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the present Invention. 

[0032) Fig. 5 Is a schematic, exploded view of a 3D diagnostic Image 90 of heart 24, following registration of the 3D 
image with geometrical map 38, In accordance with a preferred embodiment of the present Invention. TOs view Is 
generated at a bullseye rendition step 56 In the method of Rg. 2. The bullseye rendition of Image 80 comprises a stack 
5 of parallel slices 92, which are perpendicular to axis 72. The slices are preferably taken at a fixed slice Increment one 
from another along the axis. Each slice shows a section 94 of image 90, at a distance R from apex 74 that Is determined 
by the slice number. 

[0033J Fig. 6 shows slices 92 of Image 90 arrayed side-by-side, Illustrating extraction of diagnostic data from the 
slices for application to map 38, in accordance with a preferred embodiment of the present Invention. Referring, for 
10 example, to slice number 5, sectional image 94 comprises three essential parts: an inner region 100, showing the 
inside of the ventricle; a wall region 1 02, showing the myocardium; and an outer region 104, external to the heart. The 
diagnostic information of interest is in region 102. Assuming image 90 to be a SPECT image, showing perfusion in the 
heart wall, for example, region 102 will typically have the highest value of perfusion. 

[0034] At a coloration transfer step 58, the diagnostic Information from each slice 92 is transferred to map 38. Each 

is slice has a known value ot distance R from apex 74. For each angle a within the slice, point 82 on surface 80 of the 
map (Fig. 4) Is assumed to be the point at that angle that Is located radially In the middle of region 102. In the case 
that Image 90 Is a perfusion Image, point 82 Is simply taken to be the point of highest perfusion at the given angle. In 
other Imaging modalities, finding region 102 Is, for the most part, similarly straightforward. The value of the diagnostic 
data at each point 82 is preferably represented as a color applied to the corresponding region of map 38. 

20 [0035] Fig. 7 is a schematic representation of a colored geometrical map 11 0, as produced at step 58, In accordance 
with a preferred embodiment of the present invention. Because of the limited ability of a line drawing to convey qualities 
of a color image, only two different color regions appear on map 110: a well-perfused region 112, and an ischemic 
region 114. Preferably, the ischemic region has a darker or "cooler* color than the well-perfused region. In actual 
applications, in which display 36 comprises a color monitor, a broad range of different colors Is used in map 110 lo 

25 describe different levels of perfusion or of other diagnostic qualities. 

[0036] Preferably, system 20 is operated by user 22 to carry out an invasive therapeutic procedure, guided by map 
1 1 0, at a therapeutic step 60. In the present example, laser console 49 is operated to irradiate ischemic region 1 1 4 via 
catheter 30 with high-intensity laser radiation, as described in the above-mentioned PCT patent application PCT/ 
IL97/00011. The laser creates revascularization channetsin the myocardium, which are marked by system 20 with 

30 spots 116 on map 110. The combination of the imaging, mapping and therapeutic modalities enables the user to con- 
centrate the treatment in the region of heart 24 that Is known to need it, and to ensure that the region is fully covered. 
Other local therapeutic and diagnostic procedures can similarly benefit from the guidance provided by map 110. 
[0037] Although preferred embodiments are described hereinabove with reference to heart 24, the principles of the 
present Invention may similarly be applied to imaging, mapping and treatment of other organs and body structures. It 

3$ will thus be appreciated that the preferred embodiments described above are cited by way of example, and that the 
present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of 
the present invention includes both combinations and subcombinations of the various features described hereinabove, 
as well as variations and modifications thereof which would occur to persons skilled In the art upon reading the foregoing 
description and which are not disclosed In the prior art. 

40 

APPENDIX A 

[0038] This appendix provides details of step 54 in the method of Fig. 2, in which two 3D representations, P and O, 
of a chamber of heart 24 are brought into mutual registration. To begin, a rough estimate is found for the transformation 
<5 between Pand O, either manually, or using the principle axis of a bounding ellipsoid or by principle component de- 
composition. The bounding ellipsoid technique is further described In the above-mentioned European patent application 
EP 0 974 936 

[0039] Fine registration between P and O Is then preferably found using a variation of the Iterative Closest Point 
(ICP) algorithm. This algorithm is described by Besl and McKay in "A Method for Registration of 3D Shapes,' published 
50 in iEEE Transactions on Pattern Analysis and Machine Intelligence A 4(2):239-256 (1992), which is incorporated herein 
by reference. The following step3 ere repeated until convergence: 

1 . Nearest point search: For each point p in P find the closest point q on O. One can take a subset of points of P 
to improve computation speed. Simllarty, alt points ol Q can be covered, too, to ensure robustness. 
55 2. Compute registration: Evaluate a transformation T that minimizes the sum of squared distances between pairs 

of closest points (p,q). The transformation is preferably either rigid, similarity, affine or projective as described below. 
3. Transform: Apply the transformation T to all points in P. 
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(O04OJ Given two surfaces, P and O, and two sets of points, {p, G P} . fa, e 0} , step 2 of this algorithm seeks 
a transformation, T, from a family of transformations (according to the possible famine's* escribed below) that minimizes 
the mean square error, e, between the corresponding sets: 



J =1 2>-r( Pi f 



i = l 



Art Ine and Projective Transformations 

[004>.]c^ For affine transformations, defined as T[p) +Ap+t t 



15 
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A = 



a 21 *21 a 23 
*32 *33>/ 



is a 3x3 matrix, and t = 



t 2 



is a translation vector. 
We must minimize 
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we have three systems of equations: 
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7 = 1,2,3 



[0042) Let a singular value decomposition of X be X = UD\A. It then follows that a, = VOL/ty 

[0043) Projective transformations are evaluated in a similar way to the affine case. 
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Similarity and Rigid Transformation 

[0044) Whereas in affine transforms triangles are transformed to triangles, similarity transformations preserve pro- 
portions. We seek a scaling factor, c, a 3x3 rotation matrix, ft and a 3-dimenslonal translation vector, t. such that T 
5 (p) = cRp + f, wherein the error 



1 

^{R, T, c) = - £ \a - (cR Pi + t\ 



10 1 = 1 



is minimized. 

(0045) A suitable method for finding the desired similarity transform is described by Umeyama, in "Least-Squares 
Estimation of Transformation Parameters Between Two Point Patterns/ published in IEEE Transactions on Pattern 
is Analysis and Machine Intelligence, 13(4): 376-380 (1991), which Is incorporated herein by reference. Define the center 
ol mass of both P and O. 



1 n 



1 

[0048] Then define the variance of the points on both Pand O t : 



[0047] The covariance matrix between the two surfaces is 

T 



i = 1 

50 [0046] Let a singular value decomposition of E w be = l/DV^ and 



5 = 



I if det(rj)det(v) = 1 

diag(U-l) if det(rj) det(v) = -1 
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[0049] The rotation, translation and scaling of the transformation are then given by: 

R = USV 7 

t = A q ~ cR// p 

c = -— tracefos) 
a p 

to 

wherein the trace of a matrix is the sum of its diagonal elements. 

[0050) In the case of rigid transformation no scaling is applied, so that c = 1 . 



is Claims 

1 . A method for mapping a structure in a body of a subject, comprising: 

capturing a three-dimensional (3D) image of the structure comprising diagnostic information; 
20 generating a 3D geometrical map of the structure using a probe inserted into the structure; 

registering the image with the map, such that each of a plurality of image points in the image is identified with 
a corresponding map point in the map; and 

displaying the map, such that the diagnostic Information associated with each of the image points Is displayed 
at the corresponding map point. 

25 

2. A method according to claim 1 , wherein the diagnostic information: 

is related to blood flow in the structure; 
comprises local perfusion data; 
30 comprises metabolic data; 

Is related to uptake of a substance in tissue of the structure; or 
Is related to motion of the structure. 

3. The method of claim 1 or claim 2 wherein generating the geometrical map comprises bringing the probe Into contact 
35 with the structure at a multiplicity of locations on the structure, and recording position coordinates of the probe at 

the locations. 

4. The method of claim 3, wherein recording the position coordinates comprises determining the coordinates using 
a position sensor in the probe. 

40 

5. The method of any one of claims 1 to 4, wherein registering the image with the map comprises applying a trans- 
formation to at least one of the image and the map so that, following the transformation, the image and the map 
have a common axis and a common scale. 

<3 6. The method of claim 5, wherein registering the Image with the map comprises dividing the Image Into a plurality 
of parallel planar slices, perpendicular to the axis and mutually spaced along the axis, and wherein the plurality of 
image points are located In the slices. 

7. The method of claim 6, wherein registering the image with the map comprises finding an axial coordinate of each 
so of the slices and an angular coordinate of each of the image points located In each of the slices, and Identifying 

each of the image points with the map point having the same axial and angular coordinates. 

8. The method of claim 7, wherein the structure comprises a wall defining a cavity, and wherein identifying each of 
the Image points with the map point comprises finding, at the axial and the angular coordinate, the image point 

55 that is within a section of the wall. 

9. The method of any one of claims 1 to 8, wherein displaying the map comprises coloring the map to reflect the 
diagnostic information. 
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10. The method ot any one of claims 1 to 9, wherein the structure comprises a heart of the subject, and wherein 
generating the geometrical map comprises mapping an endocardial surface in a ventricle of the heart 

11. Apparatus for mapping a structure in a body of a subject, comprising: 

an imaging device, adapted to capture a three-dimensional (3D) image of the structure comprising diagnostic 
Information; 

a probe, adapted to be inserted Into the structure, so as to generate a 3D geometrical map of the structure; 
a processor, coupled to the probe and the imaging device, and adapted to register the image with the map, 
to such that each of a plurality of image points In the image Is identified with a corresponding map point In the 

map; and 

a display, coupled to be driven by the processor to display the map, such that the diagnostic information 
associated with each of the image points is displayed at the corresponding map point. 

is 12. Apparatus of claim 11 , wherein the diagnostic information: 

Is related to blood flow In the structure; 
comprises local perfusion data; 
comprises metabolic data; 
20 is related to uptake of a substance in tissue of the structure; or 

is related to motion of the structure. 

13. Apparatus of claim 11 or claim 12, wherein, to generate the geometrical map, the probe Is adapted to be brought 
into contact with the structure at a multiplicity of locations on the structure, and the processor is adapted to record 

.'5 position coordinates of the probe at the locations. 

14. Apparatus of claim 13, wherein the probe comprises a position sensor for use in determining the position coordi- 
nates. 

jo 1 5. Apparatus of any one of claims 1 1 to 1 3. wherein the processor Is adapted to register the image with the map by 
applying a transformation to at least one of the image and the map so that following the transformation, the image 
and the map have a common axis and a common scale. 

16. Apparatus of claim 15, wherein the processor Is further adapted to divide the image Into a plurality of parallel planar 
J> slices, perpendicular to the axis and mutually spaced along the axis, wherein the plurality of image points are 

located in the slices. 

17. Apparatus of claim 16, wherein the processor is adapted to find an axial coordinate of each of the slices and an 
angular coordinate of each of the image points located in each of the slices, and to identify each of the Image 

40 points with the map point having the same axial and angular coordinates. 

18. Apparatus of claim 17, wherein, when the structure comprises a wail defining a cavity, the processor is adapted 
to identify each of the image points with the map point by finding, at the axial and the angular coordinate, the Image 
point that Is within a section of the wall. 

■o 

19. Apparatus of any one of claims 11 to 16, which Is adapted to provide a map colored to reflect the diagnostic 
information. 

20. Apparatus of any one of claims 1 1 to 1 9, wherein the structure comprises a heart of the subject, and wherein 
50 geometrical map comprises a map of an endocardial surface in a ventricle of the heart. 

21. Apparatus of any one of claims 11 to 20, and comprising a medical instrument adapted to perform a medical 
procedure on the structure guided by the diagnostic information displayed on the map. 

55 22. Apparatus of claim 21 , wherein the medical instrument is contained in the probe, which is adapted to be used to 
perform the procedure locally at locations selected on the geometrical map. 

23. Apparatus of daim 22. wherein the processor Is adapted to mark on the geometrical map the locations at which 
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the procedure was performed. 

24. Apparatus of any one of daims 21 to 23 for carrying out therapeutic procedure. 

25. Apparatus of claim 24, wherein the diagnostic Information relates to local blood flow in the structure, for carrying 
out a therapeutic procedure for improving the local blood flow. 

26. Apparatus of any one of claims 21 to 23 for carrying out a diagnostic procedure. 



w 



is 



20 



25 



30 



35 



40 



55 



11 



BNSOCCIO: <£P 1 18261 8A2J_> 



EP 1 182 619 A2 




5 



EP 1 182 619 A2 



FIG. 2 
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TRANSFER IMAGE COLORATION TO 3-D MAP 
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